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Dit  rapport  is  het  eerste  van  een  serie  van  3  rapporten  over  diffractierekening  met  behulp  van 
hoogfrequent  technieken.  Het  werk  waarover  in  deze  rapporten  verslag  wordt  gedaan  is  verricht 
door  L.J.  V.  Ewijk  tijdens  zijn  detachering  bij  het  Defense  Research  Establishment  Ottawa 
(DREO)  in  Canada  van  September  1992  tot  en  met  augustus  1993.  Het  werk  heeft  daar  onder 
andere  geleid  tot  de  ontwikkeling  van  drie  codes  die  gebruikt  kunnen  worden  op  het  gebied  van 
RCS  berekening. 

Dit  rapport  geeft  een  afleiding  van  een  veelbelovende  diffractietechniek  uit  de  literatuur.  De 
reden  dat  deze  afleiding  hier  gegeven  wordt,  terwijl  deze  ook  in  de  literatuur  te  vinden  is,  is  dat 
de  afleiding  in  de  literatuur  verspreid  is  over  verscheidene  artikelen  en  hier  integraal  gegeven 
wordt.  Bovendien  worden  hier  enkele  singulariteiten  nader  beschouwd,  en  verholpen,  terwijl  dat 
in  de  literatuur  niet  het  geval  is.  Zonder  een  volledig  inzicht  in  de  afleiding  is  deze 
singulariteitenanalyse  niet  eenvoudig  te  volgen. 

De  totale  afleiding  is  weergegeven  in  een  bijiage.  De  hoofdtekst  van  het  rapport  bevat  een 
inzichtelijk  overzicht  van  de  methode,  een  vergelijking  met  gemeten  resultaten  en  de  conclusies 
die  uit  deze  vergelijkingen  kunnen  worden  getrokken. 

In  het  algemeen  kan  worden  gesteld  dat  de  gekozen  methode  goed  voldoet  en  de  gebruikelijke 
reflectieberekeningen  sterk  kan  verbeteren.  Dit  is  vooral  zo  in  het  geval  van  verstrooiing  aan 
scherpe  randen  waarbij  de  invalsrichting  van  het  veld  loodrecht  staat  op  de  rand.  Bij  schuine 
inval  is  de  invloed  van  het  berekende  diffractieveld  minder  sterk. 

De  combinatie  van  hoogfrequent  rekenmeti  oden  die  in  dit  rapport  wordt  behandeld  is  in  staat 
om  voor  vele  situaties  dc  radardoorsnede  van  objecten  te  berekenen.  Voor  die  gevallen  waar, 
volgens  de  vergelijkingen  met  meetresuitaten  die  in  dit  rapport  gepresenteetd  worden,  de 
methode  nog  tekort  schiet  wordt  aangegeven  wat  een  mogelijke  verbetering  zou  kunnen  zijn. 
Deze  bestaat  uit  het  simuleren  van  de  "surface  travelling  wave"  door  het  meermalen  toepassen 
van  de  diffractiemethode.  De  implementatie  hiervan  wordt  in  de  overige  twee  rapporten,  zoals 
hierboven  vermeld,  behandeld. 

Het  in  dit  rapport  beschreven  werk  heeft  geleid  tot  beter  inzicht  in  de  fysische  achtergronden  van 
het  verstrooiingsproces  aan  scherpe  randen,  de  verbetering  van  een  van  de  meest  veelzijdig 
toepasbare  diffractietechnieken  en  de  implementatie  van  deze  methode  in  een  rekenprogramma. 
Met  dit  rekenprogramma  is  het  nu  mogelijk  om  de  diffractie  aan  scherpe  randen  te  berekenen 
voor  bijna  alle  invals-  en  observatierichtingen.  Slechts  66n  richting  vertoont  nog  een  singulier 
gedrag,  hetgeen  problemen  opievert  als  men  het  verstrooide  veld  wil  berekenen  in  de  richting 
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van  de  voortzetting  van  de  invallende  straal  langs  een  van  de  oppervlakken  van  de  wig  indien  die 
invallende  straal  van  buiten  de  wig  komt. 

De  toegevoegde  waarde  van  de  berekening  van  diffractieverschijnselen  ten  opzichte  van 
uitsluitend  reflectieverschijnselen  komt  het  best  tot  uiting  bij  de  analyse  van  objecten  met  lage 
radardoorsnede.  De  hierboven  genoemde  uitbreiding  zal  deze  toegevoegde  waarde  nog  doen 
toenemen. 
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1  INTRODUCTION 


The  Method  of  Equivalent  Currents  (MEC)  finds  its  origin  in  the  fact  that  any  finite  current 
distribution  yields  a  finite  result  for  the  scattered  field  when  that  contribution  is  summed  in  the 
radiation  integral  [1].  So,  when  the  proper  current  distribution  can  be  found,  the  problems  with 
axial  caustics,  as  in  the  Geometrical  Theory  of  Diffraction  (GTD),  can  be  avoided.  Furthermore, 
by  using  equivalent  currents  it  is  possible  to  determine  the  scattered  field  outside  the  directions 
defined  by  the  Keller  cone. 

The  approach  in  the  MEC  currents  is  to  assume  the  existence  of  filamentry  electric  and  magnetic 
surface  currents  in  the  neighbourhood  of  surface  discontinuities,  such  as  edges,  and  to  add  these 
currents  to  the  far-field  radiation  integral.  It  is  emphasised  here  that  the  equivalent  currents  are 
not  to  be  confused  with  real  physical  currents.  They  are  not  really  present  at  the  edge,  but  are 
merely  a  mathematical  aid  to  correctly  represent  the  scattered  field.  This  is  clear  when  we  notice 
for  instance  that  these  currents  are  dependent  on  the  observation  direction,  which  is  impossible 
for  physical  currents. 

An  overview  of  the  MEC  can  be  found  in  [IJ.  The  method  is  used  by  many  authors,  see  for 
instance  Millar  (2),  [3]  and  [4],  Ryan  and  Peters  (5],  Knott  and  Senior  [6]  and  Sikta  [7].  In  these 
articles  the  MEC  is  used  for  various  scattering  problems  and  in  some  of  them  the  results  are 
compared  to  the  results  of  other  techniques. 

Knott  and  Senior  [6]  adjusted  the  equivalent  current  prescription  by  requiring  that  a  stationary 
phase  evaluation  of  the  radiation  integral  yields  the  GTD  result  for  directions  on  the  Keller  cone. 
In  order  to  extend  the  expressions  to  general  bistatic  conditions  they  replaced  some  terms  to 
invoke  dependency  on  the  incidence  and  ob.servation  direction,  a  method  later  criticised  by 
Michaeli  (8). 

Sikta  [7]  used  the  equivalent  currents  for  flat  plates  which  he  divided  in  strips  with  width  that 
approached  zero. 

Michaeli  [8]  and  [9]  proposed  a  more  rigorous  way  to  derive  expressions  for  the  equivalent 
currents.  He  stated  that  the  usual  expressions  were  based  on  GTD-like  expressions  and  were 
therefore  only  valid  for  directions  on  the  Keller  cone.  He  then  derived  new  expressions  based  on 
the  identification  of  the  MEC  line  integral  with  the  asymptotic  edge  contribution  to  the  Physical 
Theory  of  Diffraction  (PTD)  surface  integral.  From  this  resulted  the  equivalent  currents 
expressed  in  terms  of  the  PTD  surface  currents.  These  integrals  were  then  evaluated,  resulting  in 
the  equivalent  currents  that  could  be  used  to  determine  the  scattered  field  by  an  arbitrary  edge 
[8].  Later  he  improved  his  expressions  by  using  a  different  integration  variable  in  the  evaluation 
of  the  afore  mentioned  integrals  [9).  He  also  split  his  equivalent  currents  into  fringe  equivalent 
currents,  accounting  for  the  fringe  scattered  field,  and  Physical  Optics  (PO)  equivalent  currents, 
accounting  for  the  PO  scattered  field.  The  fringe  currents,  or  non-uniform  currents  as  they  were 
named  originally  by  Ufimtsev,  are  defined  as  the  difference  between  the  total  currents  on  a 
surface  and  the  PO  currents  [10].  The  step  to  split  the  equivalent  currents  into  fringe  equivalent 
currents  and  PO  equivalent  currents  was  initiated  by  Knott  [1  Ij.  He  had  compared  the  equivalent 
currents  with  the  method  of  Incremental  Length  Diffraction  Coefficients  (ILDC),  derived  by 
Mitzner  in  1974.  Knott  pointed  out  that  the  difference  between  MEC  and  ILDC  was  caused  by 
the  fact  that  Michaeli  used  the  total  surface  currents  and  Mitzner  the  non-uniform  currents  only, 
so  their  results  differ  by  terms  due  to  the  PO  currents.  Although  Mitzner  derived  his  theory  in 
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1974,  his  report  was  only  available  in  limited  distribution  for  a  long  time,  so  Michaeli  and 
Mitzner  reached  their  results  independently. 

Because  the  equivalent  currents  are  now  known  as  fringe  equivalent  currents  and  PO  equivalent 
currents  they  are  very  suitable  to  use  in  a  high-frequency  Radar  Cross  Section  (RCS)  prediction 
program.  We  can  for  instance  always  apply  the  PO  integral  to  an  object  of  which  we  want  to 
determine  the  RCS  but  only  apply  the  method  of  equivalent  currents  to  those  edges  that  are 
expected  to  give  a  significant  contribution  to  the  total  scattering.  This  last  step  can  decrease  the 
computational  time  considerably  because  a  calculation  using  equivalent  currents  is  far  more 
computation  intensive  than  using  PO.  It  also  implies,  however,  that  a  priori  knowledge  is 
necessary  to  decide  whether  an  edge  gives  a  significant  contribution  or  not  and  such  knowledge 
is  not  always  available. 

In  chapter  2  the  method  of  equivalent  currents  is  rederived  following  Michaeli  [8],  [9]  and 
casted  in  a  form  comparable  to  other  techniques  and  with  the  same  variables  as  used  in  other 
reports.  Chapter  3  will  show  some  results  obtained  with  the  computer  program  in  which  the 
MEC  is  implemented  compared  to  measurement  results. 

The  work  described  in  this  report  has  been  done  while  the  author  was  at  the  Defence  Research 
Establishment  Ottawa  (DREO)  in  the  framework  of  the  Canadian  Defence  Research  Fellowship 
Program.  The  author  had  his  office  at  the  David  Florida  Laboratory  of  the  Canadian  Space 
Agency,  where  much  of  this  kind  of  work  is  done  for  DREO,  and  was  guided  by  Dr.  S.  Mishra. 
The  measurements  that  are  used  for  the  comparisons  in  chapter  3  are  done  at  the  DDARLING 
facility  at  DFL  by  C.  Larose. 
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2  THE  METHOD  OF  EQUIVALENT  CURRENTS 


As  mentioned  in  the  introduction  this  chapter  contains  a  derivation  of  the  equivalent  currents, 
following  the  way  Michaeli  did. 

This  is  done  by  identifying  the  non-uniform  currents,  known  from  the  PTD,  with  a  magnetic  and 
electric  current  flowing  at  the  edge  under  consideration.  An  expression  can  then  be  found  for 
these  currents  after  which  they  will  be  used  in  the  radiation  integral  in  order  to  compute  the 
scattered  field.  This  will  only  be  the  edge-diffracted  field  because  the  uniform  part  of  the 
induced  total  currents  is  omitted  explicitly  from  these  expressions.  When  one  wants  to  compute 
the  total  scattered  field  one  has  to  add  the  contribution  to  this  field  by  the  uniform  currents,  for 
instance  by  means  of  a  PO  computation. 

2.1  The  equivalent  currents 

In  the  following  we  use  the  co-ordinates  and  angles  as  shown  in  figure  2.1.  The  suffixes  to  the 
angles  are  1  and  2,  instead  of  i  and  s,  becau.se  they  are  related  to  one  of  the  planes  of  the  wedge. 
The  suffixes  i  and  s  ususally  denote  angles  related  to  the  co-ordinate  axes.  When  the  edge  is 
aligned  with  the  z-axis  these  angles  are  equal,  but  generally  they  are  not. 

When  the  suffixes  1  and  2  are  used  in  conjunction  with  any  other  parameter  they  denote  the 
facenumber  that  the  parameter  is  related  to. 


Fig  2. 1:  The  general  configuration  of  the  wedge  used  to  determine  the  equivalent  currents.  The 

angles  shown  in  this  picture  are  related  to  one  of  the  planes  that  compose  the  wedge. 

The  faces  are  numbered  I  and  2. 

When  wc  assume  electric  and  magnetic  currents  flowing  along  an  edge  C  we  can  determine  the 
scattered  field  caused  by  these  currents  from  the  radiation  integral,  see  for  instance  [13]  or  [14]: 
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V  eJc 


W(r',r)dl’ 


(1) 


where  the  vector  I  denotes  an  equivalent  current,  as  opposed  to  J.  a  physical  current.  In  this 
equation  only  integration  along  the  edge  is  necessary  because  this  is  the  only  source  for  the 
currents.  This  equation  can  also  be  expressed  as: 


^(f'.r)dr 


(2) 


where  the  unit  vector  t  denotes  the  direction  of  the  edge  under  consideration  and  Ig  and  are 
now  complex  numbers  instead  of  vectors.  The  equivalent  currents  arc  always  directed  along  the 
edge. 

We  also  know  that  the  scattered  field  can  be  expressed  in  terms  of  the  physical  currents  as: 


=  J„k,x(k^xt)V|f(r'.r)dS  (3) 

for  a  peticcilj  conductii.g  wedge  The  ::umm;>iinn  in  eq  O'*  is  taken  over  both  faces  of  the 
wedge. 

Eq.  (3)  can  be  reduced,  in  the  high  frequency  limit,  asymptotically  to  a  sum  of  ray  field 
contributions  from  isolated  stationary  phase  points  and  a  boundary  contribution  expressed  by  a 
line  integral  along  C.  This  latter  integral  should  then  be  identified  with  eq  (2). 

We  will  assume  that  the  edge  diffraction  is  caused  by  surface  currents  on  narrow  strips,  Sj  and 
S2,  along  C.  Let  O]  and  oi  be  the  vectors  from  C,  perpendicular  to  C  in  the  direction  of  the 
grazing  diffracted  ray  on  face  1  and  2  respectively,  see  figure  2.2. 
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Direction  of  iacidence 


Fig  2.2:  The  direction  of  the  gniziitg  diffracted  ray  is  denoted  by  o  j  or  o  7,  the  area  where  the 

equivalent  currents  flow  is  denoted  by  S  j  or  Ss  (Sote  that  Sj  actually  is  not  visible).  The 
suffix  ( I  or  2)  to  these  ixiraineters.  indicate  the  corresponding  plane  of  the  wedge,  see 
fig  2.1. 

We  then  write  this  integral  in  such  a  way  that  only  the  asymptotic  endpoint  contribution  for  the 
a-integration  is  taken  into  account,  i.e.  the  contribution  at  o=0,  so  the  edge.  Gieen's  function  can 
be  approximated  by  this  function  evaluated  at  the  edge  multiplied  by  exp(ikOj  k^).  This  gives  for 
the  edge  diffracted  part  of  the  scattered  field  from  eq  (3); 

sin(0.1da„dl  (4) 

where  the  summation  is  taken  over  the  two  faces  of  the  wedge. 

By  comparing  eq  (2)  and  (4)  and  some  vector  calculus  we  can  find  expressions  for  the 
equivalent  currents  in  terms  of  the  integral  of  the  current  on  the  wedge.  These  currents  are  taken 
to  be  the  fringe  currents,  or  non-uniform  currents,  only.  The  solution  is  then  found  by  an 
asymptotic  end  point  evaluation  of  these  integrals,  by  using  the  solution  of  the  canonical  wedge 
problem  and  theorems  from  complex  function  theory.  This  evaluation  is  done  in  Appendix  A. 

The  result  of  this  is  given  in  eq(5)  and  (6). 
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Nksin(0|  )sin(0, ) 

Nsin((pOU(7C-(p, )  ^  Nsin(N7i-(p,)U((l-N)7i  +  9|)  ^ 


cos((p,  )  +  |.i, 
n  - 


COS(N7I-(P|  )  +  |I, 


^in((p, ) 


71  -a. 


sin(N7t - (p,  ) 


,  ,<p,  .sin(a.)  ^^7i-a,  ^  ^cp,  ^  sin(a,) 

COS( - ^)-COS(  -  ;  •  COS( - ^)  +  COs(  — )  - 

N  N  ■  N  N 


2iJ-E, 


ksin’(0,)  cos( (p I  )  +  n I 


.sin((P|  )l'(  71  -(p, )  sin(N7i  -  (p,  )U((1  -  N)7i  +  (p, ) 


cos(N7r-(p|)  +  ^_, 


<Pi 

1  1 

- - L - + - -  ) 

N  7t-«.  (p,  n-a,  tp, 

C0,S(-  -  -  )-C()S(  -)  cos( - -)  +  cos(  — ) 

.  N  N  N  N  J 

2iH,^  ^cos(0,)  cos((p,  )U{7i-tp, )  cos(N7t  -  9,  )U((  1  -  N)7r +  (p, )  ^ 

ksin(0|  )  sin(0|)  cos(  cp,  )  +  (.i,  cos( Nti  -  cp, )  +  n. 


TT  -  a. 


sin(a, 


N  N 


Hi  cos((p,) 


7C-a, 


sinla, ) 


7i  - 


cos!"-  )  +  eos(2lL)L‘‘'"V^. 


H,  cos(Nk-(p,) 


tan(0,)  tan(0,) 

N  ' 


where  p  is  given  in  eq  (A,  11 )  as: 


_  sin(0, )  sin(0, )  cosCcp, )  +  cos(0, )  co.s(02 )  -  cos^  (0, ) 

^  srn^(0,) 


(7) 
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and  a=arccos()j).  These  equations  look  very  complex  and  cumbersome  to  use,  but  they  only 
involve  simple  goniometric  functions  that  cause  no  problems  in  a  computer  program.  The  arccos 
that  must  be  calculated  can  sometimes  result  in  a  complex  number,  however,  so  care  must  be 
taken  to  use  the  correct  expression  for  this,  see  Appendix  C. 

2.2  The  diffracted  field 

With  eq's  (2),  (5)  and  (6)  we  can  compute  the  edge  scattered  field,  or  diffracted  field,  caused  by 
the  fringe  currents.  The  radiation  integral  is  easily  solved,  because  the  only  dependency  of  Ig 
and  Ipn  on  r’  is  through  the  incident  field.  This  term  can  be  kept  under  the  integral,  while  the 
other  terms  can  be  taken  outside  the  integral.  When  writing  the  solution  in  terms  of  diffraction 
coefficients  [12]  we  get: 


'’iE  "’iff 

I  = - - D.+-^= — 

kZsin'(0,)  ‘  ksin(8, ) 


2iZH, 


k  sin(0|  )sin(0, ) 


12„ 


(S) 


with  Dg,  and  as  following  froiti  eq'.s  15)  ;ind  (6).  By  using  the  radiation  integral  and 
expressing  the  incident  field  as  we  did  in  the  derivation  of  Physical  Optics,  we  find: 


E  = 


ike 


4tcR 


2iE 


Oi/ 


D,. 


0  ,£  kZsin'(9|)  "  ksin(0,) 

2iZH 

- L),„(k  xt))  I 

c 


DJ(k%x(k^xt))  + 


D,„(k^xt)}j^e-'"''-'''cir 


(9) 


ksin(0|)sin(0,) 

The  integral  can  be  solved  easily  by  using  a  parameter  representation  of  the  edge: 


dl'=Aads,  0<s<I 
f'=  Uq  +  s  Ail 


(10) 


with  aQ  the  start  point  of  the  edge  and  Aa  the  vector  along  the  edge  with  length  equal  to  the  edge. 
After  evaluation  we  find: 


I  =  Aa  sine 


k(k,-kj- 


Ail 


(11) 


During  the  evaluation  of  eq  (1 1)  we  had  to  restrict  ourselves  to  an  observation  direction  away 
from  the  Keller  cone.  This  restriction  is  only  mathematical,  however,  and  the  sine  function  in  eq 
(11)  ensures  the  correct  result  for  these  directions  of  observation  as  well. 
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Substitution  of  eq  (11)  in  (9)  and  by  omitting  some  terms  which  would  cancel  anyway  when 
calculating  the  RCS  results  in: 


e„=4»ii 


2^^  ksin'(6|) 

2iZH 


^iZH  -  -  - 

ksin(0| ) 

D„.(k.xOi 


Otr 


ksin(0,  )sin(62 ) 

A  —  *  *  AJ 

-  *  Aa  -ik(k,-k,)(a„+ — ) 

Aasinc[k(k, -k,)— ]c  - 


(12) 


2.3  Singularities  in  the  diffraction  coefficients 

From  eq  (5),  (6)  and  (8)  the  diffraction  coefficients  can  be  found  easily.  Also,  some  singularities 
in  these  expressions  can  be  pointed  out.  This  is  done  in  Appendix  B  where  the  way  to 
circumvent  these  singularities  is  de.scribcd  tis  well.  One  of  the  mentioned  singularities  will  be 
shown  here  as  well,  as  an  example. 

From  the  expression  for  D^,  we  can  find  that  this  cociiicient  becomes  singular  when: 


cos((p,)  =  -H, 

cos(Nk-(P,)  =  -(.i; 


cos(^^ — ^)  =  cos(— ) 

N  N 

K-a,  cp, 

cos( - =-)  =  -COS(--!-) 

N  N 


(13) 


When  the  observation  direction  is  on  the  Keller  cone  [  1  ],  these  singular  directions  coincide  with 
the  shadow  and  reflection  boundaries  of  face  1  and  face  2  of  the  wedge  respectively,  see 
Appendix  B.  When  the  observation  direction  is  outside  the  Keller  cone  the  singularities  occur  as 
well,  however,  so  it  is  not  possible  to  change  to  different  diffraction  coefficients,  which  would 
be  the  easiest  solution. 

The  way  to  prevent  these  singularities  is  by  evaluating  when  m  approaches  the  value  given 
by  eq  (13)  infinitely  close.  This  evaluation  is  given  in  Appendix  B  for  all  singularities,  given  by 
eq  (13),  as  well  as  for  all  the  singularities  in  Dg  and  D^.  The  resulting  expressions  are  used  in 
the  developed  computer  program  whenever  the  combination  of  incidence  and  observation 
direction  is  chosen  such  that  the  coefficient  would  become  singular. 

After  the  removal  of  the.se  singularities  only  one  singularity  will  remain.  This  is  called  the 
Ufimtsev  singularity  by  Michaeli  [8]  and  it  only  occurs  when  computing  forward  scattering. 
This  singularity  will  therefore  only  rarely  pose  problems  when  calculating  RCS. 
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3  COMPARISON  OF  COMPUTED  AND  MEASURED  RESULTS 


With  the  computer  program  that  applies  the  MEC  as  derived  in  chapter  2  it  is  possible  to 
compute  the  diffracted  field  of  general  edges.  These  edges  must  be  defined  in  a  suitable  way  of 
course,  but  there  is  no  restriction  on  the  external  wedge  angle  of  the  edge.  Any  value  between  0 
and  2n  can  be  used.  The  diffracted  field  can  be  computed  as  function  of  aspect  angle,  frequency 
or  bistatic  angle.  The  edge  is  assumed  to  be  perfectly  conducting. 

To  show  the  capabilities  of  the  method  a  set  of  objects  is  chosen  for  which  measured  data  was 
made  available.  The  improvement  by  the  MEC  on  the  scattering  computation  by  PO  will  be 
made  clear  as  will  the  shortcomings  of  a  single  diffraction  theory  in  general. 

The  objects  used  for  this  analysis  are: 

•  square  plate 

•  square  rod 

•  equilateral  triangular  plate 

•  pyramid 

•  circular  cylinder 

The  size  of  the  objects  is  such  that  they  can  be  treated  by  high-frequency  methods  in  the 
frequency  range  that  is  used  in  the  measurement  facility  at  DEL. 

The  graphs  shown  in  this  chapter  display  RCS  values  as  function  of  the  incident  angle,  usually 
the  azimuth  angle.  The  vertical  axis  displays  the  RCS  in  dB  m^,  the  horizontal  axis  displays  the 
incident  angle.  The  scale  that  is  cho.sen  to  be  consistent  for  each  object  under  consideration  so 
the  various  results  for  any  object  can  be  compared  easily.  The  pictures  always  contain  3  graphs, 
PO  results  only,  PO+MEC  results  and  measured  data  respectively. 

The  co-ordinate  system  as  used  in  the  computations  is  shown  in  fig  3.1.  When  the  orientation  of 
the  objects  is  discussed  in  the  next  sections,  the  axes  as  shown  in  this  figure  will  be  referred  to. 


Fig  3. 1:  The  co-ordinate  system  as  used  for  the  computations  in  this  chapter.  The  angles  6  and 

(p  denote  the  spherical  co-ordinates  and  are  used  as  elevation  angle  and  azimuth  angle 
respectively.  In  case  of  vertical  polarisation,  the  electric  field  vector  is  directed  as  Ey.  in 
case  of  horizontal  polarisation  as  £/,. 

In  fig  3.1  the  x-,  y-  and  z-axis  are  given  as  a  right-hand  co-ordinate  system.  With  horizontal 
polarisation  is  meant  that  the  electric  field  vector  is  perpendicular  to  the  direction  of  incidence 
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and  the  z-axis,  in  fig  3,1.  Vertical  polarisation  indicates  that  the  electric  field  vector  is 
directed  perpendicular  to  the  direction  of  incidence  and  Ef,,  as  shown  by  Ey  in  fig  3.1. 

3.1  Square  plate 

The  square  plate  used  in  this  section  has  a  sidelength  of  17.78  cm  (7")  and  measurements  have 
been  performed  using  various  frequencies,  from  5  to  18  GHz  and  for  both  linear  polarisations. 
Results  for  5  GHz  and  10  GHz  will  be  shown,  also  for  both  linear  polarisations,  displaying  the 
capability  of  the  high-frequency  methods  to  compute  results  over  a  wide  frequency  range.  For 
the  lowest  frequency  used  the  plate  only  is  about  by  3X,  just  within  the  high-frequency 
region.  The  reason  that  no  measured  results  are  given  at  a  higher  frequency  than  10  GHz  is  that 
the  plate  is  not  completely  in  the  far  field  anymore  at  those  frequencies,  which  would  distort  the 
comparison. 

During  measurements  and  computations  the  plate  is  located  in  the  yz-plane  with  its  vertical 
edges  in  the  direction  of  the  z-axis. 

RCS  (dB  m2) 


Fig  3.2:  Measured  and  computed  results  for  a  square  plate.  Frequency  used  is  5  GHz. 

polarisation  is  V-V.  The  variable  parameter  in  this  case  is  the  azimuth  angle. 
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RCS  (dB  in2) 


Fig  3.3:  Measured  and  computed  results  for  the  square  plate.  Frequency  used  is  5  GHz, 

polarisation  is  H-H. 

When  comparing  fig  3.2  and  fig  3.3  it  immediately  becomes  clear  that  the  computed  results  for 
PO  and  MEC  don't  show  any  dependency  on  polarisation.  This  phenomenon  is  known  for  PO  for 
all  objects.  The  MEC  only  displays  this  kind  of  behaviour  when  applied  to  a  flat  plate,  in  all 
other  ca.ses  different  results  wilt  be  found  for  the  two  linear  polarisations.  The  measured  data 
shows  a  broad  and  relatively  high  lobe  at  near  grazing  incident  angles  for  horizontal 
polarisation.  This  lobe  is  caused  by  a  travelling  wave,  see  [1],  a  high-frequency  scattering 
phenomenon  that  cannot  be  modelled  accurately  by  any  single  diffraction  theory. 


RCS  (dB  m2) 


Fig  3.4: 


Measured  and  computed  results  for  the  square  plate.  Frequency  used  is  10  GHz, 
polarisation  is  V-V. 
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Fig  3.4  shows  a  comparison  of  a  computation  with  PO  and  MFC  with  measured  data  at  a 
frequency  of  10  GHz.  The  size  of  the  plate,  6  X  for  this  frequency,  assures  us  to  be  well  in  the 
high-frequency  region.  Therefore,  the  computed  results  compare  better  to  the  measured  data  than 
they  did  in  fig  3.2.  The  computational  result  of  PO+MEC  differs  only  a  few  dB  from  the 
measured  data  at  maximum. 

RCS  (dB  id2) 


Fig  3.5:  Measured  and  computed  results  for  the  square  plate.  Frequency  used  is  10  GHz. 

polarisation  is  H-H. 

Fig  3.5  shows  a  similar  deviation  in  the  results  of  the  computations,  when  compared  to  the 
measured  result,  at  low  grazing  angles  as  fig  3.3.  Again  this  difference  is  caused  by  the 
travelling  wave.  Despite  this  rather  large  defect  the  improvement  of  applying  PO+MEC  instead 
of  PO  alone  is  obvious. 

3.2  Square  rod 

In  this  paragraph  the  measured  and  computed  results  for  a  square  rod  are  given.  Comparisons 
between  results  for  5  GHz  and  15  GHz  and  for  both  linear  polarisations  will  be  given.  Again, 
like  in  the  case  of  the  square  plate,  all  dimensions  are  chosen  similar  for  all  pictures  so  as  to  ease 
the  comparison.  The  size  of  the  square  rod  is  7.1  cm  by  7.1  cm  by  17.78  cm  (2.8"  x  2.8"  x  7"). 
Although  it  is  not  clear  from  the  presented  results  it  is  noted  that  the  measurements  are 
performed  with  the  rod  lying  down  whereas  the  computations  are  performed  with  the  rod 
standing  upright.  A  difference  that  only  results  in  the  two  linear  polarisations  being  changed 
from  measurement  to  computation.  The  plots  show  the  results  as  function  of  elevation  angle 
(note  that  this  is  related  to  the  computations),  so  at  0°  incident  angle  the  top  of  the  rod  is 
illuminated  perpendicularly. 
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RCS  (dB  102) 


Fig.  3.6:  Measured  and  computed  results  for  the  square  rod.  Frequency  used  is  5  GHz, 

polarisation  is  V-V. 

In  figure  3.6  we  can  see  that,  although  the  top  of  the  rod  only  measures  1  "k  squared,  the  high- 
frequency  techniques  predict  the  RCS  at  these  incident  angles  remarkably  well.  At  higher 
incident  angles  some  deviation  is  noticeable,  up  to  5  dB  at  a  peak  level,  but  just  as  we  saw 
before,  the  MEC  improves  the  PO  prediction  at  those  incident  angles  where  PO  shows  large 
discrepancies. 


RCS  (dB  m2) 


Fig  3.7 


Measured  and  computed  results  for  the  square  rod.  Frequency  used  is  5  GHz, 
polarisation  is  H-H. 


TNO- report 


FEL-94-B195 


Page 

18 


In  the  case  of  horizontal  polarisation,  fig  3.7,  the  usual  lobe  structure,  that  is  found  in  almost  all 
RCS  plots  of  simple  shapes,  is  absent  in  the  measured  data.  The  difference  between  the  PO  and 
PO+MEC  result  shows  that  the  MEC  takes  care  of  the  removal  of  this  lobe  structure  that  is,  of 
course,  present  in  the  PO  result.  It  can  easily  be  understood  that,  especially  in  the  low  incident 
angle  region,  MEC  has  a  large  influence  on  the  total  result,  because  a  leading  edge  is  illuminated 
by  an  electromagnetic  field  with  its  electric  filed  vector  directed  along  the  edge.  This  shows  in 
fig  3.7  as  well. 

RCS  (dB  m2) 


Fig  3.8:  Measured  and  computed  re.sulis  for  the  square  rod.  Frequency  used  is  15  GHz. 

polarisation  is  V-l'. 


0  >0  20  30  40  50  60  70  10  90 

elevalioo  ugle  (deg) 

Fig  3.9:  Measured  and  computed  results  for  the  square  rod.  Frequency  used  is  15  GHz, 

polaruation  is  H-H. 
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Figures  3.8  and  3.9  are  similar  comparisons  as  3.6  and  3.7.  The  frequency  used  now  is  15  GHz, 
so  we  can  be  sure  the  use  of  high-frequency  methods  is  valid.  This  is  clear  as  well  from  the 
results.  MEC  improves  on  PO  as  it  did  before.  The  eye-catching  difference  in  lobing  structure 
between  fig  3.8  and  3.9  is  again  caused  by  the  surface  travelling  wave,  as  explained  in  section 
3.1.  Note  that  the  MEC,  although  it  doesn't  take  into  account  this  surface  travelling  wave, 
predicts  the  RCS  reasonably  well. 

3.3  Equilatenil  triangular  plate 

The  equilateral  triangular  plate  in  this  section  has  a  side  length  of  17.78  cm  (7")  and  is  located  in 
the  yz-plane  during  measurements  and  computations.  Results  are  presented  for  5  GHz  and  15 
GHz  for  both  linear  polarisations. 

RCS  (dB  m2) 


Fig  3. 10:  Measured  and  computed  results  for  the  equilateral  triangular  plate.  Frequency  used  is  5 

GHz.  polarisation  is  V-V 
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RCS  (dB  in2) 


Fig  3. 13:  Measured  and  computed  results  for  the  equilateral  triangular  plate.  Frequency  used  is 

15  GHz,  palari.sation  is  H-H. 

For  fhe  lower  as  well  as  the  higher  frequency  the  MEC  doesn’t  influence  the  total  result  as  much 
as  it  did  in  the  results  presented  in  the  sections  above.  A  reason  for  this  might  be  the  fact  that  the 
edges  at  which  the  field  diffracts  are  slanted,  as  opposed  to  upright  as  they  were  before,  which 
decreases  the  magnitude  of  the  diffracted  field.  The  results  for  H-H  polarisation,  fig  3.11  and 
3. 13,  emphasise  the  need  for  a  method  to  compute  the  surface  travelling  wave. 

3.4  Pyramid 

The  pyramid  that  is  considered  in  this  section  has  a  ba.se  length  of  10  cm  and  the  angle  between 
the  faces  and  the  bottom  is  40°.  This  makes  the  height  approximately  4  cm.  Because 
measurement  results  are  only  available  at  10  GHz,  this  means  that  the  pyramid  does  not  fulfil  the 
high-frequency  requirements. 

The  bottom  of  the  pyramid  was  placed  in  the  xy-plane  with  the  tip  pointing  towards  the  positive 
z-axis. 
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RCS  (dB  m2) 


Fig  3. 14:  Measured  and  computed  results  for  the  pyramid.  Frequency  used  is  10  GHz. 

polarisation  is  V-  V. 


RCS  (dB  m2) 


Fig  3. 15:  Measured  and  computed  re.sults  for  the  pyramid.  Frequency  used  is  10  GHz, 

polarisation  is  H-H. 

From  the  results  presented  in  fig  3.14  and  3.15  it  will  be  clear  that  the  methods  used  in  this 
report  are  completely  unable  to  predict  the  .scattering  from  this  pyramid.  The  reason  for  this  is 
not  obvious,  but  there  are  a  few  possibilities.  Because  the  peak  levels  at  0°,  45°  and  90°  are 
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predicted  quite  vveli,  it  is  not  likely  that  one  of  the  methods  used  fails  for  this  object.  It  is  more 
probable,  however,  that  different  phenomena,  like  the  surface  travelling  wave,  for  instance,  play 
an  important  role  in  the  case  of  scattering  by  a  pyramid.  Another  possibility  might  be  the  fact 
that  the  size  of  the  pyramid  is  not  sufficient  to  allow  high-frequency  approximations  to  be  used. 
It  cannot  be  stated  with  certainty  what  the  reason  for  the  discrepancy  is,  but  in  my  personal 
opinion  the  latter  reason  is  the  more  obvious  one. 

3.5  Circular  cylinder 

The  cylinder  used  has  a  length  of  17.178  cm  (7")  and  a  radius  of  7.1  cm  (2.8").  It  is  located 
along  the  z-axis  with  its  centre  at  the  origin  of  the  co-ordinate  axes.  Like  in  the  case  of  the 
square  rod,  see  section  3,2,  the  cylinder  was  located  differently  during  the  measurements,  but 
again  this  does  not  pose  a  noticeable  difference. 


RCS  (dB  m2) 


Fig  3. 16:  Measured  and  computed  results  for  the  circular  cylinder.  Frequency  used  is  5  GHz. 

polarisation  is  V-  V. 


TNO-report 


TNO-report 


FEL-94-B195  I'-.ge 

25 


RCS  (dB  m2) 


Fig  3.19:  Measured  and  computed  results  for  the  circular  cylinder.  Frequency  used  is  15  GHz. 

polarisation  is  H-H. 

The  circular  cylinder  displays  a  siniilar  behaviour  as  does  ihe  square  rod.  This  is  not  surprising 
because  the  dimensions  are  closely  related.  The  MEC  improves  the  PO  prediction  and  even  for 
the  lower  frequency  the  total  result  agrees  well  with  the  measured  data,  although  the  results  for 
the  square  rod  at  the  lower  frequency  are  better.  Al.so  the  more  pronounced  lobing  structure  for 
the  15  GHz  V-V  results  are  predicted  to  a  reasonable  degree  of  accuracy. 

3.6  Summary  of  the  comparisons 

It  is  not  a  trivial  or  simple  conclusion  that  can  be  drawn  from  the  large  amount  of  comparisons 
presented  in  this  chapter.  For  the  square  plate,  for  vertical  polarisation,  the  improvement  of  PO 
by  MEC  is  very  good,  whereas  for  horizontal  polarisation  it  is  rather  poor.  For  other  shapes  the 
difference  between  the  two  polarisations  is  not  as  pronounced  and  in  general  the  results  of 
measurement  and  computation  are  in  agreement  to  a  comfortable  level. 

Bad  extremities  are  displayed  by  the  equilateral  plate  and  the  pyramid.  Possible  reasons  for  this 
are  given  in  sections  3.3  and  3.4. 

The  main  concern  at  this  moment  is  the  fact  that,  for  simple  shaper,  the  surface  travelling  wave 
can  give  an  important  contribution  to  the  total  scattering.  It  is  strongly  believed  that  an 
implementation  of  this  scattering  phenomenon  will  improve  most  of  the  results  given  in  this 
chapter. 

Finally  a  word  on  the  choice  of  the  objects  used  for  the  comparisons.  This  selection  is  made  out 
of  a  much  larger  set.  The  rea.son  that  these  specific  objects  are  chosen  is  that  they  display  the 
widest  possible  variety  of  shapes,  like  straight  ,  slanted  and  sharp  edges  as  well  as  flat  and 
curved  surfaces.  This  way  a  fair  and  objective  judgement  of  the  method  can  be  made,  see  section 
4.1. 
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4  CONCLUSIONS  AND  RECOMMENDATIONS 


4.1  Conclusions 

The  MEC  was  derived  in  the  same  rigorous  way  Michaeli  did.  The  expressions  found  reflect  the 
non-uniform  currents  only  and  therefore  the  computed  scattered  field  is  only  the  diffracted  field. 
The  advantage  of  this,  over  an  expression  reflecting  the  total  currents,  is  the  capability  to 
investigate  the  contribution  of  the  diffracted  field  to  the  total  scattered  field. 

The  method  is  made  almost  completely  uniform  in  that  the  last  remaining  singularities,  at  the 
shadow  and  reflection  boundaries,  are  removed  as  well.  These  singularities  prevented  the 
method  to  be  used  at  these  aspect  angles  and  were  not  treated  in  literature  before.  The  only 
remaining  singularity  is  the  .so-called  Ufimtsev  singularity,  occurring  when  one  wants  to 
compute  the  diffracted  field  in  the  direction  of  the  grazing  diffracted  ray  and  the  incident  field 
comes  from  outside  the  wedge.  This  situation  only  occurs  in  forward  scattering. 

From  the  results  shown  in  chapter  3  it  can  be  concluded  that  the  method  is  capable  of  computing 
the  diffracted  field  by  any  edge  for  all  directions  of  incidence  and  ob.servation,  excluding  the  one 
mentioned  above,  and  that  such  a  computation  can  provide  a  valuable  supplement  to  for  instance 
PO.  The  method  .seems  to  provide  better  re.sults  for  fields  incident  perpendicularly  to  an  edge 
than  otherwise,  which  is  something  that  needs  more  re.search. 

Furthermore  it  is  shown  that  the  combination  of  PO  and  MEC,  although  physically  only  valid  in 
the  high-frequency  region,  provide  valid  re.sults  for  simple  objects  as  small  as  only  a  few 
wavelengths.  The  polarisation  dependent  .scattering  is  not  modelled  very  well  by  the  MEC, 
although  better  than  by  PO  alone,  which  completely  ignores  this  dependency.  A  method  to 
improve  on  this  is  mentioned  in  section  4.2. 

4.2  Recommcndation.s 

Because  of  the  problems  that  the  MEC  encounters  at  grazing  incidence  for  parallel  polarisation 
it  is  strongly  recommended  that  the  method  will  be  adjusted  or  extended  in  order  to  overcome 
these  failures.  This  can  be  done  by  adding  the  capability  to  compute  multiple  diffraction  which 
accounts  for  the  interaction  of  the  diffracted  field  with  another  edge.  In  general  this  is  not 
possible  by  applying  the  MEC  again  at  the  other  edge,  because  the  diffracted  fields  do  not 
always  display  a  ray-optical  behaviour  at  the  other  edge  which  is  a  necessary  requirement  for  the 
application  of  the  MEC.  Other  diffraction  theories  have  the  same  problems  and  in  the  case  of  the 
UTD  this  was  overcome  by  a  spectral  analysis  of  the  diffracted  field.  Although  this  hasn't  been 
done  with  the  MEC  yet,  it  .seems  feasible  to  perform  such  a  spectral  analysis  with  the  MEC  as 
well. 

This  interaction  of  the  diffracted  fields  with  other  edges  is  a  simulation  of  the  surface  travelling 
wave,  the  high-frequency  phenomenon  that  cau.ses  the  discrepancies  of  any  diffraction  theory  at 
low  grazing  angles. 

When  such  an  improvement  is  implemented  the  combination  of  PO  and  MEC  can  be  expected  to 
be  able  to  give  accurate  results  for  all  kinds  of  objects  in  a  large  frequency  region.  Especially  the 
RCS  prediction  of  not  to  complex,  low-RCS  objects  will  benefit  from  this. 
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APPENDIX  A:  DERIVATION  OF  THE  EQUIVALENT  CURRENTS 

By  comparing  eq  (2)  and  (4)  in  the  main  text  we  see  that; 


I,(k,x(k^xt))-^I,„(txkJ  =k,x[k^xjJ^e‘'‘‘'"'‘'’dx„] 


By  taking  the  scalar  product  of  eq  ( A.  I )  with  k^xfk^xt)  or  (k^xt)  we  arrive  at  an  expression  for  Ig 
or  Ifp  respectively: 


I,  =  .  ^  k^ -j^CtxkJxJ J(asin(9,),gcos(9|))e’''*°^'*daj 

1  '  r  ^ 

=  — -T;:^J-t[k.xfj(osin(0,).ocos(e|))e‘‘^‘°‘^'»do] 
sin‘(0  J  V  e  J 


So  the  evaluation  of  Ig  and  I,,.,  can  now  be  completed  by  an  asymptotic  end-point  evaluation  of 
the  current  integral  in  eq's  (A. 2)  for  a  canonical  wedge  problem.  When  we  use  the  following 
property  of  the  3-dimensional  solution  of  the  canonical  wedge  problem,  see  James  [13): 


J(x,z)  = 


(A.3) 


and  evaluate  the  vector  products  in  eq  (A. 2)  we  arrive  at: 


.  f rj_  j.COSCP,)  J 

sin(0,)  J 


l^=-Z 


J  COS(03  )C( 
sin(0 

e' 

sin(02)-^  * 


jj^gik(ak.-cov(o.))  sin(0^)d^j 


The  currents,  J,  in  these  expressions  are  supposed  to  be  the  fringe  currents,  or  non-uniform 
currents,  only  and  not  the  total  surface  currents  and  are  evaluated  at  (asin(0|),  0). 

We  can,  in  the  usual  way,  express  the  currents  as  the  cross  product  of  the  surface  normal,  (0,  1, 
0),  and  the  magnetic  field  at  the  surface.  Then,  with  Maxwell's  equations  and  the  expressions  for 
the  incident  fields  at  the  surface: 
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VxE  =  -iCO|.lH 
VxH  =  icoeE 


(A.5) 


and 


■ik(k,  F  ) 


P  ^-ik(xsin(Q,  )cosl»?,)+2cos(0,)) 


rV  -»k(.xs»njO,)coM<p,)+zcos(6,)) 


(A.6) 


we  can  easily  find  an  expression  for  as  function  of  E^  and  which  is  needed  in  the 
expression  for  the  current.  Tliis  equation  is; 


H.  = 


I 


iksin"(0|) 


e  aE, 

~^  +  cos(e,)-^ 
'  U  dy  dx 


(A.7) 


From  James  we  now  use  the  canonical  wedge  solution  to  relate  E,  and  H7  to  Ej^  and  as: 


*12' 


E,  =  E„,,[ii(X.V}/-(p,)-ii(X,\t/  +  (p,)] 
Hz  =  Ho„[u(X,v-(p,)  +  Li(X,\|/  +  (p,)] 


with 


X  =  kasin‘(0,) 


,  sin(— )e 

dTriN^r 


^  ^^iXcoNl^) 


C  <|) 

cos(— )  -cos(— ) 

N  N 


and  <t)  the  polar  angle  of  the  field  point  in  the  xy  plane,  so  for  our  purposes  (1)=0. 

By  substituting  eq  (A.7),  (A. 8)  and  the  expression  of  the  surface  current  in  terms  of  the  magnetic 
field  in  (A.4)  we  get  the  expressions  for  the  electric  and  magnetic  equivalent  currents.  This  will 
be  done  first  for  the  magnetic  current. 

From  the  definition  of  u  following  eq  (A. 8)  we  see  that  u(X,-<t))=u(X,(t))  because  (|)  is  only  present 
in  the  cos  term.  Also,  with  the  definition  of  X,  we  can  change  to  X  as  the  variable  of  integration. 
This  yields: 


_  ^sin((p;)f 
sin(02) 


H,,  Jii(  X, -tp, )  +  u(  X ,  tp, )]  sin(0,  »da 


(A.9) 
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I 


0) 


2ZHo„sin((p,) 

ksin(0,)sin(02) 


ju(X,(p,)e-^’^ 


with  (j  following  from  the  comparison  of  eq  (A.9)  and  (A.  10)  as; 


(A.IO) 


ka(ak^  -cos'(0|)) 

^  kasin-(0,) 

,  (A.ll) 

_  sin(0|)sin(03  )cos((p,)  +  co.s(6,  )cos(6,)-cos~(6|) 

sin’(0| ) 

The  derivation  of  Ig  is  a  little  more  complex,  but  we  can  make  use  of  the  expression  for  from 
eq  (A.IO). 

In  terms  of  the  surface  currents  vve  have: 


_cos(0,)_cos(y,)  r J 
sin(0, )  •’ 


(A.12) 


The  second  term,  including  the  minus  sign,  is  equal  to  cos((p2)/(Ztan(02))  1^,  with  as  given 
in  eq  (A.IO).  We  therefore  only  have  to  elaborate  the  first  term  in  eq  (A. 1 2).  We  know  that 
can  be  expressed  in  terms  of  and  with  eq  (A.7)  and  (A. 8)  in  terms  of  Ej^  and  Hj^.  We  then 
find  for  the  z-component  of  the  surface  current: 


I  J  ,  fg  p  3u(X.-(p,)  au(X,(p,), 

^  iksin'(0,)'  ^4  ay  ay  ' 


+  cos(0,)Ho,J 


au(x.-(p,)  au(x,(p,). 


(A.  13) 


Because  X=kosin2(0j)  and  x=osin(0|)  we  can  replace  dx  by  dX/(ksin(0i)  and  also,  from 
geometrical  point  of  view: 


au  _  1  au 
ay  p  atp, 

_  ksin(0,)  au 
X  ^ 
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From  the  definition  of  u,  following  eq  fA.8)  we  can  deduct: 
8u  _  3u 

All  this  substituted  in  eq  (A.  13)  gives: 


1  ,  ksin(e,)  3u(X,(p|) 

^  X  atp, 


+  2H„^ksin(e,)cos(e,) 


au(x,(p,) 


(A.14) 


So  we  have  for  the  first  term  of  eq  f  A.  1 2): 


2i  -E„ 


I  = — IS -  __ 

ksin'  (6| )  X  0(p 


fli^e'''"dX  +  ^ 

J  Y  t- 


2iH„i,cos(e,)  f  3u 


ksin"  (0, )  •'  dX 


(A.15) 


with  p  as  defined  in  eq  (A.  1 1 ). 

The  second  term  of  eq  (A. 12)  follows  easily  from  eq  (A. 10).  see  text  following  eq  (A. 12)  and 
reads: 


Ztan((p,)  ksin(0|)sin(03) 

: ,  2H,_^cos(.(.=_)cos(9,  )  ,  ^(x.,p,)e“"'dX 

ksin(0|)sin(0,)  •' 


(A.  16) 


We  now  have  found  the  expressions  for  the  electric  and  magnetic  equivalent  currents  in  eq 
(A.15)+(A.16)  and  (A. 10)  respectively. 

The  integrals  in  these  equations  will  be  evaluated  by  using  the  method  of  stationary  phase, 
resulting  in  their  asymptotic  end-point  value  at  the  edge  [8].  The  function  u  is  defined  as  a 
contour  integral  in  the  complex  plane,  see  figure  A.l  for  the  contour. 

From  eq  (A. 10),  (A.15),  (A. 16)  and  the  definition  of  u  following  eq  (A8)  we  see  that  we  must 
compute  three  different  contour  integrals.  We  will  do  this  starting  with  the  one  that  is  given  in 
eq  (A.  10).  When  we  substitute  the  definition  of  u  in  this  equation  and  omit  the  term  in  front  of 
the  integral  for  now,  we  have: 


TNO-reporl 


(A.  19) 
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COS(^,)  =  -|I 

=  7t-arccos(|i) 


=  n-a 


The  residues  in  the  poles  given  by  eq  (A.  19)  can  be  calculated  with: 


Res= 


-2sin(^i^)sin(^^) 

2N  2N 


(A.20) 


{A.21) 


which  follows  from  elementary  goniometric  relations  and  the  solution  of  eq  (A.  19).  The  residue 
at  the  pole  given  by  eq  (A.20)  is: 


Res=  lim(q-q,) 


isin(— ) 

_ _ 

[cos{— )  -co.s(~)](-2sin( - ~r— )sin(  — — )) 


(A.22) 


-isin( - ) 

^ _ Ln _ 

[cos(^^^ — ^)cos(— )]sin(a) 

N  N 

According  to  the  residue  theorem  the  integral  evaluates  as  27ti  times  the  sum  of  the  residues 
when  the  direction  of  the  contour  is  counterclockwise  around  the  poles,  which  is  the  case  here. 
So  the  final  result  for  this  integral  (eq  (A. IS)  can  now  be  written  as: 


1  N  ^  N  ^ _ I 

N  N 


(A.23) 


The  remaining  two  integrals  are  computed  in  the  same  manner. 

The  differentiation  of  u  to  X  in  eq  (A.  15)  can  be  performed  without  difficulty  and  results  in  an 
equation  that  resembles  eq  (A.  17)  very  much.  This  equation  can  be  treated  exactly  as  is  done 
above  and  results  in. 
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L  = 


2iE„ 


+ 


NkZsin'(e,) 

cos(0,) 


Nsin((p,)U(7t-(p,) 


sin(— ) 
N 


C0S((P,)  +  |I, 


cos(- 


L)-COS(^) 

N  N  , 


Nksin‘(0.) 


sin( - ~) 

N 


COS((p,)U(Tt-(p,) 

cos((p,)  +  ^,  cos(^^)-cos(^) 

N  N 


2iH„„cos((p,; 


Nksin(0,)  tan(0,; 


sin( - 

N 


1 


NU(n-(Pi) 

cos«p, )  +  n,  ■"  )  Sin(a, ) 

N  N 


(A.28) 

Eq  (A. 27)  and  (A.28)  give  the  equivalent  currents  flowing  on  face  1  of  the  edge.  The  currents  on 
face  2  can  be  found  by  performing  the  following  substitutions  in  these  equations; 


(p,  =  NK-(f), 

0,  =n-0,  (A.29) 

z  =  -z 

We  can  then  determine  the  equivalent  currents  on  the  second  face  of  the  wedge  and  subtract  this 
from  the  value  of  the  first  face  in  order  to  find  the  complete  expressions  for  the  equivalent 
currents.  These  expressions  are  given  in  the  main  text,  in  eq's  .(5)  and  (6). 
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APPENDIX  B:  REMOVING  SINGULARITIES  FROM  THE  EQUIVALENT 

CURRExNTS 


In  this  chapter  we  will  show  how  the  diffraction  coefficients,  derived  in  the  main  text,  can  be 
made  uniform.  The  way  they  are  given  now  still  contains  non-uniformities  that  become  apparent 
at,  for  instance,  the  optical  boundaries.  Although  this  can  be  circumvented  quite  easily  in  a 
computer  program,  it  is  far  more  appealing  to  have  expressions  that  can  be  used  for  every  angle 
of  incidence  and  observation. 

This  is  done  in  this  chapter  by  determining  the  limiting  value  of  every  expression  that  becomes 
singular.  In  the  computer  program,  where  these  coefficients  are  used,  one  must  switch  between 
the  original  expression  and  the  limiting  one  when  an  optical  boundary  is  approached.  Because 
the  original  expressions  only  become  singular  at  the  optical  boundary  itself,  the  mentioned 
switch  can  be  made  as  close  as  0.01  or  0.001  degree  from  the  boundary,  depending  on  the 
numerical  accuracy  of  the  computer  u.scd. 

First  the  expressions  for  the  diffraction  coefficients  are  given  again  and  the  singularities  are 
pointed  out.  Further  we  look  more  specifically  at  the  terms  that  actually  become  singular.  In  the 
main  text  we  can  find  Dg  and  and  slightly  rewritten  they  read: 

^  _  sin((p,)U(TC-(P))  ^  ■sin(N7i-(p,)U((l-N)7t-KP| 
cos((p,)-t-|.i|  COS(N7t-(P|)  +  |.t, 


sin(  ■ 


Tc-a 


+  ■ 


N 


sin((p;) 


N[cos( 


sin(a,) 


N 


-)-COS(^)] 
iN 


sin(- 


,  n-a.. 


-t-- 


N 


N[cos(  ^ )  +  cos(  — )  1 

N  N 


.sin(N7r-(p,) 

sinfa,) 


(B.l) 


sin((P|)U(n-(P|)  Mn(NK-(p,)U((l-N)7i  +  (P|) 

cos((p,  co.s(N7i-(p,  )-f  jij 

sin(^)  sin(^) 

_ _N_ _ lji_ _ 

N  [  cos(  —  ■  ^ -  cos(  ^ )  1  N[cos(-  —)-Fcos(^)] 
N  N  N  N 


(B.2) 
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D  -  ^ 

cosftpi  )U(Tt;-(p, )  cos(N7t:-(p,  )U((1-  N)7r  +  <P| 

tan(0,) 

1 

cos((P|  )  +  |.i,  cos(Nrt-(p,  )  +  (i, 

cos((p,  )U(n-(p, )  cos(N7i-(p,)U((l-  N)Tr  +  (p, 

tan(0, ) 

co.sfcp,  )  +  |.l, 

-a. 

N  ' 

COS(N7C-(p,  )  +  |i, 

fl,  cos((p2) 

N 

,7i-a.  (p,  ■ 

cost - )  -  cos( — ) 

N  N  _ 

7t  -  a, 

M„(  ^  ■  ) 

sinfa, ) 

tan(0,)  tan(0j) 

|i,  cosfNir-tpj) 

_J - 

N 

Tt-a,  (p, 

cos( - =- )  +  cos(  — ) 

sinfa, ) 

tan(0|)  tan(0,) 

Starting  with  the  expression  for  we  can  see  that  this  coefficient  is  singular  when: 


(B,3) 


cos((p,)  =  -U| 
cos(NTr-(P|)  = 


cos( 

cos( 


(B.4) 


N 

n-a, 

N 


N 

)  =  -cos(— ) 

N 


From  these  equations  we  can  find  the  values  of  a|  and  (Xj  for  the  singular  directions: 


a,  =  rt-tp,  vTt  +  cp, 

'  '  '  (B.5) 

a,  =  n-Nn  +  (PiVK+N7i-(P| 

The  values  of  a|  in  the  first  line  cause  the  terms  for  face  1  to  be  singular  and  the  values  of  0-2  in 
the  second  line  cause  the  terms  for  face  2  to  be  singular. 

When  the  observation  direction  is  on  the  Keller  cone  of  directions  and  we  express  tp2  in  terms  of 
aj  or  a2,  we  see  that  these  singular  directions  coincide  with  the  reflection  and  shadow  boundary 
of  face  1  and  face  2  respectively.  In  general,  however,  these  singularities  exist  outside  the  Keller 
cone  as  well,  so  we  don't  make  any  a  priori  statement  about  (P2. 

The  components  of  that  become  singular  for  aj  =  Tt±(p|  are  taken  together  and  we  then  try 
to  determine  the  value  of  this  combination  when  tt]  approaches  one  of  these  values  infinitely 
close.  Thc.se  terms  are: 
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Sin(  - 


n- a, 


N 


)sin(cp, ) 


sin((pO 


N  sin(a,  )[cos( 


Tt-a,  (p, 

- L)_co.s(^ 

N  N 


cos((p. )  +  cos(a| ) 


(B.6) 


When  we  substitute  a|  =  n-(p|+5  in  equation  (B.6)  and  take  the  limit  for  5  approaching  0  we 
arrive  at: 


t',  =  lim 

'  5-.0 


Sint 


‘P, 


N 


).sin((p, ) 


Nsin((p-6) 


<p,-6  ffi, 

co,s(  — )  -  cos(  — ) 
'  N  N 


^  +  • 


sinftp, ) 


cos((p, )  -cos((p,  -6) 


(B.7) 


When  the  .sum  of  the  two  cosines  in  the  denominator  is  expanded  we  get  terms  that  contain  5/2. 
For  later  use  we  therefore  introduce  the.se  .same  terms  in  the  numerator  as  well  in  the  following 
way: 


sin((p, ) 

-5/2  5/2^  ]■ 

N  N  sin((p, -5/2-6/2) 

sirjftp,  -5) 

-*'"<1’'  -5/2)sin(6/2) 

L  N  N  J. 

V  u  .... 

(B.8) 

We  now  expand  the  terms  in  the  numerator  and  divide  equals  terms  out  so  that  we  reach: 


t',  =  lim- 

‘  5-»0 

sin(9, ) 

5/2 

cos(-^) 

N 

cos((p,  -6/2) 

sin((P|  -5) 

5/2 

2N  sin( - ) 

2Nsin(^‘ 

2sin((n,  -6/2) 

L  N 

N 

(B.9) 

In  the  limit  we  can  simply  set  6  equal  to  0  in  almost  all  terms,  except  in  those  where  the  term 
would  become  infinite  by  doing  so.  This  happens  in  the  two  terms  with  sin(8/x)  in  the 
denominator.  For  5  very  small,  however,  we  can  replace  sin  by  its  argument  and  cos  by  1.  When 
we  do  this  we  see  that  the  singular  terms  cancel  exactly  and  we  get  for  the  limit; 
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sin((p_, ) 
sin((p, ) 


1 


2Ntan(^) 

Is 


2lan((p,) 


(B.IO) 


When  a]  =  7i+(p|+5  is  substituted  in  cq  (B.6)  and  we  follow  the  same  derivation  we  arrive  at 
exactly  the  same  expression  for  the  combination  of  the  singular  terms. 

This  expression  should  therefore  be  used  instead  of  the  original  one  when  tt]  is  near  one  of  the 
above  mentioned  values. 

When  aj  is  near  one  of  the  values  given  in  eq  (B.5)  we  must  determine  the  limit  for  two 
different  terms  and  for  ao  =  it+NK-(p|+5  or  0-2  =  Tt-Nn+tpj+S,  with  5  approaching  0  in  both 
cases.  The  singular  terms  are  in  this  case; 


.  ,  7t-a,  ,  . 

sin( - ^  )  .sin(  Nti-  tp, ) 

N 


sin(Nix  -  tp, ) 


-^iM  cos(NT:-(P|)  +  cos(a.,) 


Nsin(a,  )[cos( - -•-)  +  cos(  — )] 

N  N 


(B.ll) 


By  Substituting  one  of  the  values  of  ao.  taking  the  limit  for  5  approaching  0.  introducing  terms 
with  6/2  in  the  numerator  as  before  and  evaluating  the  resulting  expression  we  arrive  at: 


_  sin(N7t-(p,) 
'  sin(Nn-(P|) 


1  1 

N 


(B.12) 


This  expression  is  valid  for  both  values  of  ao  that  cause  these  terms  to  become  infinite. 

By  now  we  have  expressions  for  all  singular  terms  of  D,.,.,,  so  it  is  now  possible  to  compute  the 
correct  value  of  this  coefficient  for  every  angle  of  incidence  and  observation.  Following  the 
same  procedure  for  Dg  we  find  for  the  singular  terms  of  face  1 ; 


t.,  =- 


sin(^) 

N 


sin((p|) 


N[cos(^^)-cos(^)]  ‘^os(9,)  +  cos(a,) 

N  N 


(B.13) 


and  when  the  limit  is  computed  this  results  in; 
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t'5 


1 


2Ntan(^) 

N 


1 

4* - 

2tan((p, ) 


For  face  2  the  singular  terms  are; 


(B.14) 


j  ^ _ '  N  ' _ siii(N7t-(p,) 

N(cos(^^)  +  cos(?l)|  cos(N*-(p,)  +  cos(a,) 
N  N 

and  their  limits  are; 

,  I  1 

N 

As  possible  singular  terms  in  ^ve  have  for  face  1 ; 


(B.15) 


(B.16) 


sin(- 


n-a 


^9  ~ 


N 


Nsin(a. ) 


N  N 


cos(a| ) 
tan(6, ) 


cos((p,) 
tan(0, ) 


1 

+ - - 

cos(9|  )  +  cos(a| ) 


cos(9,)  ^  cos((p,) 
tan(0,)  tan(0,) 


(B.17) 


These  terms  can  be  treated  similar  as  and  Dg  above  with  only  one  small  difference.  In  this 
case  we  arrive  at; 
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Page 

B.6 


t’,  =  lim 


5/2  q>, -5/2^ 

lim— JL_  (  N  cos((p, -5)  ^  cos((p,) 

5-osin((p, -6)  2Nsin(-^'~--^^) 

.  “  '  N  N 

cos(5/2  cos((p|  - 5/2)  Y^os((p, )  cos((pj)^ 


(B.IS) 


2sin(5/2)  2sin(9,  -  5/2)  tan(6, )  tan(9,) 


Now  those  terms  with  siii(5/x)  in  the  denominator  do  not  cancel  out  as  before.  The  reason  for 
this  is  that  they  are  not  multiplied  by  the  same  factor  although  in  the  limit,  when  6  goes  to  zero, 
these  factors  are  equal.  When  co.s((pl-5)  is  expanded  and  the  goniometric  functions  are  replaced 
by  their  small  argument  values  v«nere  appropriate  we  finally  get  for  this  limit; 


1  sinftpi )  ^  1 

sin((P|)  tan(0|)  2tan((P|) 


[cos((P|)  cos((p,) 

ONtan(^)  vtanO,)  tan(e,) 
N  ) 


(B.19) 


Again,  as  we  saw  before  u  ith  and  Dg,  this  limit  is  valid  for  both  singularities. 
The  .singular  terms  for  face  2  are: 


cosfrXj)  cos(N7i-(p,) 


Ns.n(a,)  cos(^)Acos(^)  ^^"(6,) 

N  N 


_ 1 _  COs(N7l-(p,  )  ^  COSCNTl-tPj) 

cos(N71- cpi ) +  cos(aj )  tan(0,)  tan(0,) 


(B.20) 


and  following  the  same  way  as  above,  we  find  for  the  limit  of  these  terms  approaching  one  of  the 
optical  boundaries; 
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1 

sin(  Nit  -  <P| )  1 

\ 

1  1 

1 

" coslNir - (pi )  ^  cos(Nk - (pj ) 'l 

t  J  1  ““ 

sin(Nrt -  (P| ) 

tan(0| ) 

2Ntan(^)  2tan(N7r-(p,) 

In  j 

,  tan(0,)  tan(6,)  j 

(B.21) 


which  is  valid  at  both  singular  directions. 

We  have  determined  the  values  of  all  pertinent  singular  terms  in  the  coefficients  D^,  Dg  and 
Dj^.  These  expressions,  together  with  the  original  ones,  can  be  used  to  compute  the  diffraction 
coefficients  for  all  angles  of  incidence  and  observation,  with  the  exception  of  the  Ufimtsev 
singularity.  Whenever  the  observation  direction  is  chosen  such  that  one  of  the  diffraction 
coefficients  becomes  singular  one  has  to  use  the  appropriate  limit.  A  reasonable  estimate  is  to 
change  from  one  expression  to  the  other  at  0. 1  degree  from  the  singular  direction. 
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APPENDIX  C:  COMPUTING  THE  INVERSE  COSINE  FOR  ALL 

ARGUMENTS 


When  computing  the  inverse  of  the  cosine  in  the  programs  that  apply  the  MEC  it  is  possible  that 
the  result  is  a  complex  number.  Not  every  computer  can  handle  this  in  a  correct  manner,  so  in 
Older  to  prevent  problems  we  will  use  an  unambiguous  expression  for  this  function.  Such  an 
expression  can  easily  be  derived,  as  will  be  shown  in  this  chapter.  We  start  with; 


a  =  arccos(li) 

|i  =  cos(a)  (C.I) 

=  I[e“+c"“] 

In  the  usual  way  we  can  express  sin(a)  as  a  sum  of  exponential  functions  and  combination  with 
the  last  line  of  cq  (C.  1 )  results  in: 


e'“  =}a  +  i7l-(.l' 

J  I - 

a  =  -ilnjfl  +  iVl-M  J 

As  long  as  p  lies  between  -1  and  1  the  square  root  in  eq  (C.2)  can  be  evaluated  in  a 

straightforward  manner.  When  the  value  of  p  goes  beyond  these  . .  however,  one  has  to  be 

careful.  In  accordance  with  Michaeli  (8),  we  make  the  following  choice  for  a  in  those  cases; 


a  =  -iln{p  +  7p-  -l},(l<p) 

a  =  -iln{p  -7fl'-  <(.t) 

Eq  (C.2)  will  be  used  in  the  other  cases. 

With  this  choice  of  sign  we  are  assured  that  for  any  p; 


(C.3) 


sin(a)  = 


(C.4) 
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